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ABSTRACT 

Cold fronts (CFs) — density and temperature plasma discontinuities — are ubiquitous in cool cores 
of galaxy clusters, where they appear as X-ray brightness edges in the intracluster medium, nearly 
concentric with the cluster center. We analyze the thermodynamic profiles deprojected across core 
CFs found in the literature. While the pressure appears continuous across these CFs, we find that 
all of them require significant centripetal acceleration beneath the front. This is naturally explained 
by a tangential, nearly sonic bulk flow just below the CF, and a tangential shear flow involving a 
fair fraction of the plasma beneath the front. Such shear should generate near-equipartition magnetic 
fields on scales < 50 pc from the front, and could magnetize the entire core. Such fields would explain 
the apparent stability of cool-core CFs and the recently reported CF-radio minihalo association. 
Subject headings: cooling flows — galaxies: clusters: general — hydrodynamics — intergalactic 
medium — magnetic fields — X-rays: galaxies: clusters 



1. INTRODUCTION 

In the past decade, high resolution X-ray observations 
have revealed an abundance of density and temperature 
discontinuities known as cold fronts (CFs). They are 
broadly classified as merger CFs (not discussed here) and 
cool c ore CFs; for review, see [Markevit ch fc Vikhlininl 
(|200l . 

Core CFs are observed in more than half of the other- 
wise relaxed, cool core clusters (CCs, iMarkevitch et all 
I2003T) . in distances up to ~ 400 kpc from the center. 
They are usually nearly concentric or spiral, and mul- 
tiple CFs are often observed in the same cluster. The 
\ temperature contra st across such a CF is q = T Q /Ti ~ 2 
(|Owers et al.l 120091 ). where inside/outside subscripts i/o 
refer to regions closer to/farther from the cluster cen- 
ter, or equivalcntly below/above the CF. The plasma 
beneath the CF is typically denser, colder, lower in en- 
tropy and higher in metallicity than the plasma above 
it. Usually, little or no pressure discontinuity is identi- 
fied across these CFs, suggest ing negligible CF motion 
(|Markevitch fc Vikhlirlir3 l2007h . 

Such CFs were modeled as associated with large 
scale "sloshing" oscillations of the intracluster 
medium (ICM) driven , for example, by mergers 
(|Markevitch et al.l l200ll hencefort h M01), possibly in- 
volvi n g only a dark matter subha l o dTittlev fc Henriksenl 
120051: lAscasibar fc Markevitchl 120061) . or by weak 
shocks /acoustic waves displacing cold central plasma 
(jChurazov et al.ll200l iFuiita et al.ll2004[ ). 

We argue that a tangential shear flow across and be- 
neath CFs is ubiquitous in CCs. One line of argu- 
ment relies on thermodynamic profiles across CFs found 
in the literature. We show that all these profiles re- 
veal a distinct deviation from hydrostatic equilibrium, 
and argue that this can be naturally interpreted only 
as shear. An independent indication for the presence of 
shear across CC CFs is their remarkable sharpness and 
stability, which implies strong magnetic fields along the 

1 Einstein fellow 



CF plane. In the absence of radial CF motion, the most 
natural explanation for the persistence of such fields is 
shear across the CF. 

In we show that deviations from hydrostatic equi- 
librium, previously reported in two CC CFs, are com- 
mon among such CFs, and argue in $3] that this directly 
gauges bulk tangential flows and shear along and beneath 
these CFs. In $4] we show that such shear can magne- 
tize the core, and produce near-equipartition fields along 
CFs, thus stabilizing them. We summarize and discuss 
the implications of our model in $5] We assume a Hub- 
ble constant H = 70 km s _1 Mpc -1 . Error bars are la 
confidence levels. 

2. DEVIATIONS FROM HYDROSTATIC EQUILIBRIUM 

The mass density p(r) and pressure P(r) radial pro- 
files near CFs in RX J1720+26 and A1795 were reported 
to be inconsistent w ith hydrostatic equilibrium (M01; 
iMazzotta et al.l 120011 ). Indeed, assuming such equilib- 
rium, a r = —p~ 1 d r P — GMg/r 2 = 0, where a r is the 
radial acceleration and G is Newton's constant, would 
imply an unphysical gravitating mass profile M g (r) that 
becomes abruptly larger just outside the CF. In RX 
,11720, M g derived from the above assumption jumps by 
a factor of 5 (at a 2.5tr confidence level) at the CF, and 
there is marginal evidence for a pressure discontinuity 
which translates t o radial CF motion w ith Mach num- 
ber M = 0.4l°;^ (jMazzotta et al.ll200lh . In A1795, M g 
jumps by a factor of 2 (2.5cr) and the pressure appears 
continuous across the CF (Figure [TJ inset; see M01). 

To investigate the prevalence of this phenomenon, we 
calculate the pressure acceleration discontinuity 

Aa r = a r i — a r o = A( — - — ) = — — A(AT) (1) 

P fJ-r 

across various CFs for which deprojected thermodynamic 
profiles have appeared in the literature. Here we defined 
A A = Ai— A a (for any quantity ^4) and A = d In P/d In r 
(typically A < 0), with ks being the Boltzmann constant, 
p ~ 0.6m p being the average particle mass (assumed con- 
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stant), and m p being the proton mass. The gravitating 
mass M g (r) is assumed to be continuous everywhere. 

We extract all the thermodynamic profiles across CFs 
in cool cluster cores from the literature. We select only 
profiles that have been deprojected along the line of sight, 
where at least two radial bins exist on each side of the 
CF such that d r P can be estimated. When possible, we 
omit radial bins which are very far (> 4 times in radius) 
from the CF, to avoid excessive contamination. We end 
up with six CF profiles in which AT can be calculated 
based on 2-3 radial bins on each side of the front, supple- 
mented by three CFs for which only the density has been 
deprojected. We present the (dimcnsionlcss) normalized 
acceleration discontinuity 8 = {c 2 /r)~ l Aa r in Figure [TJ 
along with the literature references. Here, c = (TP/p) 1 / 2 
is the sound velocity, with T = 5/3 being the adiabatic 
index. 




r [kpc] 

Figure 1. Normalized discontinuity <5 in pressure acceleration, 
plotted agai nst radius r for C Fs in vari ous CCs (filled sym- 
bols). Data: ITanaka et at] l200l. for A496). I Johnson et al.1 120101 . 
for A1644), M01 ( f or A1 795), ISanders et all H2003 , for A2204), 
and ISanders et ai l 112009, for 2A0335) . Data with p but no 
T deprojection: Mazzotta lrtaLl 12001] . for RX J1720+26) and 
Markcvitch ct al. (20Q0, for A2142). The pressure is (consistent 
with being) continuous across all these CFs. Also shown are the 
acceleration discontinuities in some deprojected regions where no 
CF is present (empty cyan symbols), used as a control sample. The 
— 1 < 8 < (shaded) region is roughly consistent with subsonic 
shear and flow beneath the CF, < A(u 2 ) = (ui + ii )Aii < c 2 . 
Inset: Gravitating mass (in Mq) assuming hydrostatic equilibrium, 
as a function of r (in kpc) in two opposite sectors in A1795, one 
harboring the CF analyzed (adopted from M01). 



All the Aa r values we find are negative, and inconsis- 
tent with zero at confidence levels ranging between 1.1 
and 5.4a. The CF data points shown in Figure 1 arc 
consistent with a best fit S = — 0.78± 0.14, i.e. nonzero 
at 5.5cr (the average is not necessarily physically mean- 
ingful). As a sanity check, we apply the exact same 
procedure to parts of the same deprojected thermody- 
namic profiles where no CF is present. This is possible 
in a few cases, where the profiles extend sufficiently far 
from the CF, shown in Figure fT]) as empty cyan sym- 
bols. These continuum 6 values are consistent with zero 
or are slightly positive (within 2.5er from zero; a resolu- 



tion effect explained below). Overall, they are best fit by 
S = +0.07±0.08. This indicates that the a r discontinuity 
at CFs is real. 

Furthermore, the effect is qualitative and robust. As T 
jumps by a factor of q ~ 2 as one crosses outside the CF, 
hydrostatic equilibrium would require a slower radial de- 
cline of pressure outside the CF, Xi/X ~ q. However, the 
pressure profile typically steepens outside the CF. Hence, 
two phenomena contribute to the Aa r < discontinuity: 
the density drop, and the pressure steepening. Neither 
of these effects can be inverted by spherical deprojection. 
The differentiation procedure used to estimate A removes 
long wavelength errors, so the significance of S is proba- 
bly better than that estimated above. Finally, note that 
resolution effects introduce an opposite bias, toward pos- 
itive Aa r , as (— p~ l d r P) typically declines with increas- 
ing r. This accounts for the slight tendency of continuum 
regions toward Aa r > 0. 

We conclude that deviations from hydrostatic equilib- 
rium are common in CC CFs, with less pressure support 
below the CF, Aa r < 0, such that if the plasma above the 
CF is in equilibrium, the plasma below should be falling. 
This phenomenon, previously reported in RX J1720 and 
in A1795, has been interpreted as evidence that the core 
is sloshing, with plasma near the CF radially accelerat- 
ing in some phase of the oscillation cycle. For the appar- 
ently stationary CF in A1795, cool plasma would then 
be sloshing at its maximal radial displacement, where it 
has zero radial velocity but nonzero radial acceleration 
(M01). 

We note, however, that in this picture, without ad- 
ditional effects, the discontinuity in radial acceleration 
Aa r = a r: i — a r ^ < would cause an unphysical gap 
to open along the CF. A different interpretation of the 
observations is therefore required. Below we argue that 
the most natural interpretation of the data involves a 
noncontinuous (within the observational resolution) cen- 
tripetal acceleration of flow along the curved CF. 

3. TANGENTIAL SHEAR FLOW BENEATH CFS 

Consider an infinitesimal, two-dimensional piece Ti of 
the CF plane, and an inertial "CF frame" S' instanta- 
neously comoving with Ti perpendicular to the CF plane, 
as illustrated in Figure [21 Let T\ be an arc in Ti parallel 
to the local flow beneath the CF. We may approximate 
T\ as a segment of a large circle with radius R. Con- 
sider the spherical coordinate system S' = (r',8,4>) with 
T\ lying at radius r' — R and colatitude 6 = tt/2. Its 
origin, the center of the circle, is typically close but not 
necessarily coincident with the cluster center. 

The r' component of the velocity u measured in S' 
must vanish at J-2, u r > = 0, as ?' is perpendicular to the 
discontinuity. The r' component of the acceleration is 
continuous across J^, A(dtU r ') = 0. Hence, taking the 
difference in the r' component of the Euler (momentum) 
equation between the two sides of Ti yields 



A(u 2 )/r' = Aip^dr'P) 



(2) 



where the pressure P = Pth + Put may include thermal 



and nonthermal components, u 



u • u 



ui, and 



we assumed slow changes in the CF pattern, u • Vw r / <C 
u 2 /r' . Note that although u is confined to the CF plane, 
it does not have to be parallel on the two sides of T<i. 
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Figure 2. Local geometry of the CF and the inner flow. The 
radius of curvature R of the CF profile T\ (solid curve) at a given 
position (circle) is typically slightly larger than the distance r mea- 
sured along r (dot-dashed) from the cluster center. The flow (filled 
arrows, arbitrarily shown as an outflow) near the CF must parallel 
Ti in a frame S' (dashed axes) comoving with the CF. The ob- 
served M g profile implies that shear flow extends far beneath the 
CF (filled arrows of different sizes). 

Equation ([2]) indicates that the observed discontinuity 
A(p~ 1 d r 'Pth) > must be balanced by a discontinuously 
larger centripetal acceleration u 2 /r' beneath the CF. The 
discontinuity cannot be plausibly attributed to P n t, be- 
cause the continuous P t h, noncontinuous ^ Pth combina- 
tion observed would then require predominantly nonther- 
mal pressure at la rge, sometimes r > 50 kpc radii, which 
is unlikely (e.g., IChurazov et al.ll2008D . Henceforth we 
neglect P nt . 

We conclude that shear flow across CFs is the most 
natural interpretation of the observations. Thus 

u? - ul = A(u 2 ) ~ (k B /fi)A(X'T) ~ -rAa r > , (3) 

where A' = 3 In P/dlnr' ~ (r'/r)(f • r')A ~ A. Hence, 
CFs are tangential discontinuities, with shear and faster 
tangential flow beneath the CF, Ui > u a . In particular, 
Au < Au2 < Mi, where Au2 = [A(w 2 )] 1//2 is measurable. 

In terms of the local sound velocity, Equation <j3j) be- 
comes 

TA{u 2 ) = A'Jc 2 = (\> - X' q)c 2 . (4) 

Typically the P slopes are Xi > —0.7, steepening to 
A < — 1 above the CF, so a subsonic flow in the CF 
frame («; < Cj) requ ires q < 2.4. Ind eed, q < 2.5 is 
observed in CC CFs (jOwers et al.ll2009f ). Assuming the 
subsonic flow and continuous pressure, Equation (J4j) also 
imposes an upper limit on the velocity above the CF, 
Ml = [uo/cof < q- 1 + X'JT - Xi/(Tq). 

As an example, consider the CF in A1795. The M g 
discontinuity was previously interpreted as reflecting ra- 
dial acceleration, a~2x 10~ 8 cm s~ 2 . Interpreted as 
tangential flow, this corresponds to A«2 ~ (ar) 1 / 2 ~ 
750 km s _1 , uncertain to within a factor of ~ 2. This 
velocity difference constitutes a considerable fraction of 
the sound velocity, Au 2 — 0.67ci ~ 0.55c o . It places a 
lower limit A4i > 0.67 on the (tangential) Mach number 
beneath the CF, and an upper limit Am < 0.67ci on the 
(nonzero) shear across it. Assuming A4i < 1, it also im- 
poses an upper limit M. < 0.61 on the flow above the 
CF. 

Note that the tangential flow inferred from Aa r / 
may coexist with perpendicular (radial) motion of the CF 



through the ICM, inferred from an appar ent pressure dis- 
conti nuity, for example in RX J1720+26 (Mazzott a et al.l 
l200l . 

The M g (r) profile near a CF can be compared to other, 
equidistant regions in the core, preferably on the oppo- 
site side of the cluster, that do not harbor discontinu- 
ities. In both cases examined, A1795 (Figure [TJ inset) 
and RX J1720+26, M g (r > r ) is similar on both sides 
of the cluster above the CF radius ro, but inertia ap- 
pears to be missing (M g is too small) in most of the 
volume beneath the CF. In the tangential flow picture, 
the missing inertia region beneath the CF corresponds 
to an extended tangential bulk flow. Consider again the 
CF lying r ~ 85 kpc south of the center of A1795. 
Here, inertia appears to be missing in a considerable, 
0.7 < r/ro < 1 region, involving ~ 1/2 of the gas in that 
sector. There is no evidence for inertia discrepancy below 
~ 0.6rrj. Hence, plasma within the discrepant radii range 
cannot be uniformly flowing; shear must be present. 

The existence and extent of the shear flow region are 
also gauged by nearly radial optical (and sometimes X- 
ray) filaments that lie below some CFs. In A1795, a pair 
of fila ments extend out to ~ 0.6ro ([McDonald fc Veilleuxl 
120091 and references therein), just beneat h the shear re- 
gion. In M87, the South- West filaments (|Forman et al.l 
|2007[ ) extend out to ~ 0.6r . They are nearly radial out 
to <~ 0.4ro, but have some curvature at 0.4ro < r < 0.6ro 
which may be driven by shear. The (sometimes curved) 
termination of filaments just below the shear region sug- 
gests that they may have been eroded by the shear. 

4. MAGNETIC SHEAR AMPLIFICATION 

An independent argument implying shear across CFs 
stems from their remarkable sharpness and stability. Sta- 
tionary CFs arc stable against Raylcigh- Taylor insta- 
bilities, as the inside plasma has lower entropy. How- 
ever, particle diffusion, heat conduction, and instabil- 
ities such as Kclvin-Helmholtz (KHI) and Richtmyer- 
Mcshkov, could potentially broaden and deteriorate the 
CF on short, subdynamical time scales. In contrast, ob- 
servations indicate that CFs are not only stable over 
cosmological time scales, but in fact the transition in 
thermodynamical properties subtends a small fraction of 
the Coulomb mean free path of proto ns (jVikhlini n et al.l 
I2001al iMaFkevitch fc Vikhlininl 120071 ). 

This suggests that transport across CFs is magneti- 
cally suppressed. Such suppression is expected in CFs 
that move through the ICM, due to draping of mag- 
netic fields and possibly the formation of a magnetic bar- 
rier (plasma depletion layer) , as found in obser vations of 
plane tary magnetospheres and in simulations (|Lvutikovl 
120061 and references therein). However, the processes 
that suppress transport in the more common, stationary 
type of CFs, were thus far unknown. Here we argue that 
shear magnetic amplification near CFs naturally sustains 
the strong magnetic fields parallel to the CF needed to 
suppress radial transport. 

Consider first the limit where CFs are infinitely sharp 
tangential discontinuities. The magnetic field perpendic- 
ular to such a CF must vanish, and there is no shear mag- 
netic amplification. However, in this limit the stability 
against KHI inferred from obser vations requires a strong 
magnetic field B, as shown by [Vikhlini n et ail ()2001bD 
for the merger CF in A3667. Namely, the shear veloc- 
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ity differ ence must be smaller th an an effective Alfvcn 
velocity (lLandau fc Lifshitzlll96(jl §53), 



AuY B 2 3(1 + q) 2 



< 



4;rp 5/?i 



(5) 



and in addition 



(B 4 x B ) 2 > 2irp (B< x Au) 2 + (B x Au) 2 , (6) 



where B 2 = (B 2 + B 2 )/2, p = 2p lPo /{ Pl + Po ), and p = 
Pth/P§ = 8wnkBT I B 2 , with n being the particle number 
density. 

Equation (|5|) imposes a lower limit on B, 



B > B„ 



2irpi 



\Au\ 



\Au\ 



12(1 + ?) a 



B, 



50 ev cm J 



(7) 



The large shear velocities derived in §21 thus imply large 
_B, up to ~ 40% of its equipartition value B eq for which 
,8 = 1. An asymmetry Bi/B ^= 1 would require B well 
aligned with Au. In the case Bi = B = B, oti ^ a a , 



/ 2(sin Oj + sin a ) 

\l '. 2~7 \ L>min j 

y sin {on - do) 



(8) 



with a being the angle between B and Au, so stronger 
fields are needed if misaligned with the shear. 

A magnetic, plasma depleted barrier may form. As 
long as its temperature does not greatly exceed T , it re- 
mains strongly magnetized. Note that the corresponding 
density drop could possibly be observed more easily here 
than in merger-type CFs, as there is no confusion with a 
stagnation region. As shown below, equipartition fields 
are expected only very close to the CF, so P n t can still 
be neglected at > kpc distances, as assumed in $3] 

How did strong magnetic fields form along stationary 
CFs, and how are they sustained over cosmological time 
scales? Perpendicular transport suppression could pre- 
serve such fields practically indefinitely in the absence 
of large perturbations {e.g., mergers) and disruptive ki- 
netic plasma effects. Hence, they could be remnants of 
an early stage where the CF was formed. But if such 
initial fields are absent or decay, shear driven magnetic 
amplification can naturally generate or replenish B^, by 
stretching B r i =/= magnetic structures advected by the 
flow via line freezing, d(H/p)/dt ~ [(B/p) ■ V]u. 

Consider a weakly magnetized, B <C B eq stage where 
the kinematic viscosity v is a fraction /„ of its Spitzer 
value. The velocity discontinuity across such a CF is 
gradually smoothed out, forming a transition layer of 



thickness I 



(i/t) 



1/2 



7tf' 2 Trn:n*/Gyr) 1/a kpc 



beneath the CF at time t, where T4 = (Tj/4 keV) and 
?i_2 = (rii/10~ 2 cm -3 ). This enables rapid shear mag- 
netic amplification, on an e-fold time scale r ~ l/u. 
Equating the viscous and amplification times suggests 
that a magnetization layer of thickness 



Nv/u 



5QNf v MT l T'tn 



ln_\ pc 



(9) 



develops over a very short time scale 

t - N 2 v/u 2 ~ 4 x lO 4 ^ 2 /^ -2 ^ 372 ?!.:^ yr , (10) 

with N being the growth factor. Typically 7V/„ < 1; 
N ~ 1 may suffice in the case of a weakening magnetic 
barrier. 

The high B estimates above hold only in the vicinity 
of CFs, where fast shear is measured. As mentioned in 
§31 evidence suggests that the associated shear involves 
a fair fraction of the plasma beneath the CF. Shear am- 
plification may thus effectively magnetize the entire core, 
at least beneath the outermost CF. For a simple model 
where the shear is constant in a layer of thickness A 
beneath the CF, the energy of the magnetic field is am- 
plified in proportion to [B^j B r /) 2 , where 



B<t> .a 

B r > 



1/2 



A 



( 1 



10 kpc / V 108 Y r 



(11) 



Saturation of this process depends on the flow details. 
For example, in a spiral flow there is a characteristic cut- 
off time when magnetic structures leave the shear region. 

We have seen that CFs are strongly magnetized, and 
that the associated shear flow can magnetize much of 
the plasma beneath them. Magnetic processes should 
therefore be enhanced below, and in particular near 
CFs. This naturally explains the recently discov- 
ered coincidence be tween CFs and the edges o f ra- 
dio minihalos (MH) (Mazzotta &: Giacintuccil 120081 ); see 
IKeshet fc Loebl l|2010h . 



5. DISCUSSION 

We interpret deviations from hydrostatic equilibrium 
previously observed near CFs as evidence for centripetal 
acceleration. This reveals a tangential shear flow, nearly 
sonic at the CFs and extending well beneath them. Such 
behavior is shown to be common among CC CFs, as illus- 
trated in Figure [T] The measured discontinuity A (it 2 ), 
given by Eqs. ([3j|4]), constrains the shear and inside flow, 
Am < Au2 < Ui, and imposes an upper limit on M. . 
The apparent CF stability against KHI suggests near- 
cquipartition magnetic fields along the CF (Eqs. [THHl; 
this could manifest in a shear-amplified, < 50 pc mag- 
netic layer. Shear magnetic amplification could rapidly 
magnetize the entire core up to the CFs, and explain 
some of the observed properties of radio MHs. 

In our model, core CFs are part of extended tangential 
discontinuities seen in projection. They directly trace 
both the orientation and approximate magnitude of bulk 
flows in the ICM. CFs are usually nearly concentric, so 
B should be nearly tangential and the polarization ra- 
dial. Note that this magnetic configuration differs from 
the predictions of magnetic co mpression models , where 
B is typically isotropic {e. g., iGitti et aD 120021 ) or ra- 
dial (jSoker k, Sarazinlll990l) . but is similar to the satu- 
ratio n of heat flux d riven buoyancy instabilities in the 
core (lQuataertj|20"0"l . In MHs,' we predict detailed mor- 
phological correlations with X-ray gradients, which trace 
projected CFs, and radio polarization parallel to these 
gradients. 

Such fast, extended flows would modify the energy 
budget of the core and could alter its structure, for cx- 
amp le by distributing heat and relaxing the cooling prob- 
lem. iChurazov et al.l (|2003f ) pointed out that the absence 
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of resonant features in the X-ray spectrum in Perseus 
indicates motions with at least Ai = 1/2 in the core. 
Although initially interpreted as turbulent motion, this 
may reflect the bulk motions discussed here. The pres- 
ence of such flows could be tested, for example using 
high resolution spectrosc opy and future X-ray polariza- 
tion measurements (e.g., iSazonov et al.ll2002|) . 
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